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In this paper a detailed numerical study (in frames of the Slonczewski formalism) of magnetization 
oscillations driven by a spin-polarized current through a thin elliptical nanoelement is presented. We 
show that a sophisticated micromagnetic model, where a polycrystalline structure of a nanoelement 
is taken into account, can explain qualitatively all most important features of the magnetization 
oscillation spectra recently observed experimentally (S.I. Kiselev et al., Nature, 425, 380 (2003)), 
namely: existence of several equidistant spectral bands, sharp onset and abrupt disappearance of 
magnetization oscillations with increasing current, absence of the out-of-plane regime predicted by 
a macrospin model and the relation between frequencies of so called small-angle and quasichaotic 
oscillations. However, a quantitative agreement with experimental results (especially concerning the 
frequency of quasichaotic oscillations) could not be achieved in the region of reasonable parameter 
values, indicating that further model refinement is necessary for a complete understanding of the 
spin-driven magnetization precession even in this relatively simple experimental situation. 

PACS numbers: 85.75.-d, 75.75.+a, 75.40.Gb, 75.40.Mg 

I. INTRODUCTION 

Magnetic excitations and magnetization switching in nanostructures induced by a spin-polarized current (SPC), 
first predicted theoretically 1 and soon discovered experimentally are now one of the most intensively studied topics 
in solid state magnetism due to their importance both from the fundamental point of view and for numerous possible 
technical applications in, e.g., microwave generators and MRAM cells (see recent reviews 3 ). In the past decade a 
substantial progress has been achieved by the elaboration of analytical theories explaining the spin-transfer driven 
magnetization dynamics and by refining experimental techniques in order to enable quantitative experimental studies 
of these phenomena. 

Analytical models employing the scattering matrix formalism™ and taking into account spin-dependent reflection 
and transmission of electrons, precession of electron spins and the structure of Fermi surfaces for concrete materials^ 
have been developed which allow the calculations of the spin-transfer torque arising due to the current crossing the 
normal metal-ferromagnet interface. The derivation of the density operator corresponding to the normal metal spacer 
(N) in a FM/N/FM structure combined with a two-current (for spin-up and spin-down electrons) model enabled the 
evaluation of a spin torque in a complete trilayer structure^ which was successfully tested and extended using the 
Boltzmann equation formalism 9 . This way an asymmetric dependence of the spin torque on the angle between the two 
FM layers could be derived&LS whereby the parameters governing this asymmetry depend on the conductivities and 
geometry of the multilayer device under study. A large effort has also been made to establish the relation between 
th e spin- transfer phenomena, interlayer exchange coupling and Gilbert damping in multilayer structures (see Ref. 
flolllllfT^ll^ and citations therein). 

At the same time first high quality experimental results were obtained by studying magnetization oscillations 
induced by a spin-p olarized current in various systems. Magnetization dynamics in nanopillars was investigated in 
Ref. Ib3lllalllfll7lllll Detailed studies of the steady-state magnetization oscillations were reported also for the point- 
contact geometrjii2i22i2i. In these papers several different precession modes were identified and the dependencies of the 
precession frequency and microwave oscillation power on the external field and dc current strength were investigated. 

Quantitative comparison between theory and experiment remains, however, a delicate issue due to fairly compli- 
cated remagnetization processes involved. Experimentally accessible multilayered nanostructures have thicknesses in 
the range of several nanometers and the lateral sizes about 100 nm and larger which is far above the critical size 
for the single-domain behaviour. For this reason one expects the formation of relatively complex domain structures 
during the magnetization switching, especially taking into account that for most equilibrium magnetization states 
the spin-induced torque has different directions for different parts of the nanoelement (see, e.g., the discussion in 
Ref. |22|) . Indirect evidences for such a domain- mediated switching can be found already in first semiquantitative 
experimental studies^, where a quasistatic switching of a nanoelement by changing the external field or the electric 
current strength was investigated. Recent experiments have confirmed that a macrospin (single-domain) approxima- 
tion can not explain several important observations concerning the magnetization dynamics in nanopillars and point 
contactsi^i&Si, although the macrospin model was able to predict some of its qualitative aspects like the existence 
of a steady-state precession regime by itself (for the most detailed study in frames of a macrospin approach see Ref. 
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0). 

For this reason the full finite-element micromagnetic modelling is required in order to understand whether the dis- 
crepancies between theoretical predictions and experimental observations are due to a poor theoretical understanding 
of a basic physics involved or because the formation of complicated domain structure significantly alters the experi- 
mentally measurable system properties. Until the present time very few such studies have been carried out. In the 
pioneering papeipSS devoted to this subject the formation of complicated magnetization states already in a nanoele- 
ment with lateral sizes as small as 125 x 125 nm was found. In the latter detailed studj^ of the magnetization 
dynamics of a 64 x 64 nm square with the thickness 2.5 nm it was predicted that the behaviour of this element is 
nearly single-domain. However, the results of Ref. |2f| turned out to be at least partially incorrect for the parameters 
of magnetic materials typically used experimentally: it was shown^SiSi that for a nanoelement with this thickness sig- 
nificant deviations from a single-domain behaviour appear already for lateral sizes b « 30 nm and a completely chaotic 
magnetization pattern was found for b > 60 nm, which clearly emphasizes the importance of full-scale micromagnetic 
simulations. 

Up to our knowledge, systematic micromagnetic studies concerning the spin-injection driven magnetization dynam- 
ics have been performed only for a SPC induced noise in spin- valves^. In addition, in several short reports22i2ii some 
qualitative aspects of the magnetization dynamics of elliptical nanoelements (concerning mainly the transition from 
a homogeneous to a non-coherent magnetization precession) were discussed. 

In this paper we present a systematic study of the magnetization dynamics of an elliptical nanoelement under 
the influence of a spin-polarized current. The paper is organized as follows: In Sec. [H]we outline our simulation 
methodology, paying special attention to the justification of our choice of the simulation parameters. In Sec. IIIII 
we present our simulation results, starting with a 'minimal' micromagnetic model and carrying on with adding the 
effects of a polycrystalline structure and thermal fluctuations. In the same section we analyze also the influence of the 
exchange stiffness constant on the magnetization dynamics, keeping in mind that among various magnetic parameters 
of a ferromagnetic material the experimental determination of this constant is the most difficult task. In Sec. II VI we 
compare our results first with numerical simulations of other groups and then with the experimental studies of the 
microwave oscillations in such a nanoelement. At the end of this last Section we discuss how the factors which were 
not included in our model could affect the magnetization dynamics of the system under study. 

II. NUMERICAL SIMULATIONS: METHODOLOGY AND CHOICE OF MATERIAL PARAMETERS 

Numerical simulations were carried out using our package MicroMagus^i, whereby the spin injection was included in 
form of the Slonczewski torque T = (aj/Ms) ■ [M x [M x S]] (S is the spin polarization direction of the current through 
a nanoelement). The time evolution of the magnetization configuration was obtained by integrating the Landau- 
Lifshitz-Gilbert equation of motion for the system magnetization using the modified Bulirsch-Stoer algorithm. This 
method includes an adaptive step size control, what is especially important when simulating dynamics of a strongly 
non-homogeneous magnetization configuration. 

Keeping in mind the intended comparison of our simulation results with those obtained experimentally, we have 
made an effort to choose the system parameters as close as possible to their values for experimentally studied na- 
noelements. Among these experimental studies we have selected the publication of the Cornell groupie, where ones 
of the qualitatively most interesting results on Co/Cu/Co system were reported. 

An elliptical element with the geometry corresponding to that used in Ref. fl4l i.e., with lateral sizes 130 x 70 nm 
and thickness d = 3 nm was simulated. The element was discretized in plane using the mesh size 2.5 x 2.5 nm (we 
have checked that further mesh refinement did not lead to any noticeable changes in the simulation results). 

The choice of the magnetic material parameters requires special justification. 

The saturation magnetization of thin Co layers, as it is well known, may be different from that of the bulk material. 
For this reason we have used the saturation magnetization 4nMs = 10 kG as measured in Ref. ^] on a Co/Cu 
multilayer system where each Co layer was 3 nm thick (see below also our analysis concerning the influence of the 
saturation magnetization value on the final results). 

The value of the exchange stiffness constant A for a bulk Co and thin Co film is also a subject of controversial 
discussions. Some earlier measurements have lead to the values in the range A = (1 — 2) x 10 -6 erg/cm (see Ref. 
1321 and citations in Ref. l33l pp. 368-381). In latter experiments^ performed using the inelastic neutron scattering 
on bulk Co crystals and recent measurements^ carried out on thin Co films using the elastic small-angle neutron 
scattering significantly higher values A — (2.8 — 3.1) x 10~ 6 erg/cm have been found. For this reason we have chosen 
the value A = 3.0 ■ 10 -6 erg/cm as a basis and have also studied how the smaller values of A would affect the results. 

The treatment of the Gilbert damping (represented in the standard version of the LLC-equation by the double vector 
product [M x [M x H eff ]] multiplied by a constant factor A) is at present a subject of an intensive research and is also 
closely related to the problem of a correct choice of the spin torque term. First, even in the absence of an externally 
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driven spin-polarized current the value of the dissipation constant in thin film systems may substantially deviate from 
its value for the corresponding bulk material (see a short survey of experimental data collected in Ref. ITol) . This 
effect is usually attributed to a spin transfer (and subsequent relaxation) from a ferromagnet into a normal metal 
layer— and/or to a similar process leading to a "dynamic exchange coupling" between two FM layers^. In addition, 
it was showni£ that due to the dependence of the spin pumping and the spin accumulation on the angle between 
the magnetization of FM layers in a FM-N-FM systems the damping depends on the instantaneous magnetization 
configuration of FM layers. 

The same mechanisms are responsible for a non-trivial dependence of the magnitude of the spin torque T on the angle 
9 between the FM layer magnetizations (the trivial part of this dependence, V ~ sin 6, is contained in the double vector 
product [M x [M x S]]). Various theoretical approaches2i2i2& have led to the conclusion that in a symmetric spin valve 
the functional dependence of the spin torque on 9 can be cast into the form T(9) ~ PA 2 sin 9 /[(A 2 + 1) + (A 2 — 1) cos 9} , 
where P is proportional to the difference between the resistances for the electrons with opposite spin directions and 
A may be expressed as a function of a total resistance and conductance of a non-magnetic metallic spacer. For a 
non-symmetric spin valves more complicated formulae were derived^. 

Both topics discussed above - (i) dependence of the dissipation parameter A on the average and local characteristics 
of a multilayer and (ii) a complicated dependence of a spin torque magnitude on the angle between the magnetization 
vectors in adjacent layers - can and should in principle be incorporated into micromagnetic simulations. However, the 
study of the influence of these factors on the magnetization dynamics is for obvious reasons a separate and difficult 
problem. In our opinion, this problem should be addressed after the magnetization dynamics driven by a spin-polarized 
current is sufficiently well understood for the model where both the damping parameter and the coefficient before 
the Slonczewski torque term are assumed to be constant. As we shall see below, this is still not the case if one goes 
beyound the macrospin approximation, i.e., when the complicated magnetization structures arising in nanoelements 
under realistic conditions are taken into account. For this reason we have decided to put both A = Const and 
a j = Const (so that T(0) ~ sin9) in order to find out (i) which qualitative features of the experimentally observed 
steady state precession and switching of magnetic nanoelements can be explained within this simplified treatment of 
the spin-transfer induced torque and (ii) whether this model can provide quantitative agreement with experimental 
data for reasonable parameter values. In future, the comparison of our results presented here with simulations 
employing more sophisticated spin torque models will allow the unambiguous identification of dynamic phenomena 
arising due to the abovementioned non-trivial properties of the Gilbert damping and the complicated ^-dependence 
of the spin torque magnitude. 

Within the approach outlined above we are left with two problems: how to choose the value of the Gilbert damping 
A and the proportionality coefficient between the experimentally measured current values / and the Slonczewski term 
magnitude aj for our system. We are not aware of any independent measurements of A for the experimental situation 
similar to that studied in Ref. HtL which are reliable enough to be used as a simulation input. Hence we have adopted 
the following strategy to choose this parameter. First, it is well known from general considerations and confirmed by 
numerical simulations^ that the current threshold I m - m for the onset of the steady state magnetization oscillations 
growth linearly with A. At the same time, we have found that the critical current value / ma x for which the oscillations 
disappear is nearly independent on the damping. For this reason we have chosen the damping which enabled us 
to reproduce a ratio of the maximal and minimal currents J m ax/^min for which a significant microwave power was 
observed experimentally' 1 ' 4 . This procedure lead to the value A « 0.04. This damping is significantly higher than the 
usually reported values of A in Co (A ~ 0.01 — 0.02) which might indicate that the effects discussed above lead (on 
average) to a substantial enhancement of the damping in the system under study. 

Establishing of the relation between the total current I (or the current density j) and the spin torque magnitude aj 
is required for the calculation of the current induced magnetic field (Oersted field). Following the discussion above, 
we did not try to calculate the corresponding proportionality coefficient using various available theoretical models. 
Instead, we have again used for this purpose the value of the current threshold I m - ln reported in Ref. Q. Namely, 
we have calculated the Oersted field for the critical current value I m i n (which is w 2 mA for iT ex t = 2 kOe in Ref. 
Il4^ corresponding to the onset of steady state oscillations. Afterwards, keeping the Oersted field constant, we have 
increased the value of aj until the steady state precession emerged thus obtaining the critical value of the spin torque 
a c J. The ratio between I m ; n and a c J was then used as a constant proportionality coefficient between aj and I (we 
need the latter to calculate the Oersted field for the known geometry) when we increased the spin torque value aj in 
our studies. 

We have chosen the coordinate system with Oxz plane coinciding with the element plane and Ox axis parallel to 
the long ellipse axis. All simulation results presented here were obtained for the external field Ho with components 
Hq x — 2000 Oe, Hop = 20 Oe, Hq z = 100 Oe. The -ffox-value corresponds to that used for several measurement series 
reported in Ref. Il4l the small in-plane deviation H$ z was introduced to mimic a corresponding small angle between 
the long ellipse axis and the external field 14 , and the very small out-of-plane component Hq v was introduced to remove 
an artificial numerical degeneracy present when the field is applied exactly in the layer plane. The spin polarization 
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vector S was chosen to be antiparallel to the external field thus assuming that the magnetization direction of the 
polarizing (lower) layer coincides with the external field (the upper layer magnetization dynamics is driven by the 
electrons reflected from the polarizing layer so that their average spin polarization is opposite to that of the lower 
layer). 

For a direct comparison of simulation results with experimental data we need to calculate the spectral power 
of the resistance oscillations. To do this, we make an assumption that in the CPP-geometry (current perpen- 
dicular to plane) the resistance variation due to the magnetization dynamics for the current flowing through 
the i-th discretization cells of the free (magnetization direction mf ce ) and fixed (mf x ) layers can be written as 
ARi = Ai? max (l — cos#i)/2 = AR max (l — (mp ce mf x ))/2. We note in passing that at present more complicate depen- 
dencies of the magnetoresistance on the angle between the adjacent layer magnetizations are discussed (in addition 
to the papers mentioned above, see, e.g., Ref. @). However, at this stage of our research where we have to use the 
angular-independent spin torque magnitude aj, we prefer to use the standard cos-like dependence written above thus 
treating both phenomena (magnetoresistance and spin torque) in a self-consistent manner. 

For the device under consideration the variation of the resistance due to the magnetization precession is much 
lower than its value in the absence of magnetoresistive effects. Using this simplifying assumption, the total time- 
dependent resistance variation of our discretized micromagnetic model treated as a system of N c parallel resistors can 
be expressed as a function of the fixed layer magnetization direction m (this layer is assumed to be homogeneously 
magnetized) and the average magnetization direction (m frcc ) of the free layer as 



A f? 

Ai?totW = {1 ~ ' (m fr °°))} (1) 

Taking into account that the in-plane magnetization direction of the fixed layer for the geometry under study 
coincides with in-plane external field direction (given by the unit vector e h = Hq/Hq) and that the out-of-plane 
component of m is vanishingly small, we finally obtain the formula 



AiMi) = -^fH 1 - ( e "K rec > + (2) 

which enables a direct calculation of the magnetoresistance oscillation spectrum from the known time dependencies 
of the average in-plane magnetization projections of the free layer. 

For further references we also note that the frequency of the homogeneous FMR-mode of the extended thin film with 
magnetic parameters and under external conditions given above would be /^ R = ( 7 /27r) • {H {H + AttM s )) 1/2 ~ 13.7 
GHz. 



III. NUMERICAL SIMULATIONS: RESULTS AND DISCUSSION 



A. Magnetization dynamics for T — 0: a 'minimal' micromagnetic model with the reference parameter set 

We begin with the description of the steady state precession without thermal noise. In this case there exist 
for a fixed external field and increasing current (increasing aj in our formalism) a sharp transition from a stable 
magnetization state (time-independent magnetization configuration) to a steady-state regime with regular oscillations. 
In this subsection we consider a 'minimal' micromagnetic model, where the polycrystalline structure of a nanoelement 
is neglected, so that no random magnetic anisotropy of crystal grains is taken into account. 

We start the analysis of our simulation results from the regular oscillation regime, which is characterised by a nearly 
homogeneous magnetization configuration during the precession. For the nanoelement parameters (AttMs = 10 kG, 
A = 3.0 • 10~ 6 erg/cm and A = 0.04) and the external field (H = (2000, 20, 100) Oe) given above, the current 
threshold for the oscillation onset in our "minimal" model was found to be a c j = 0.3310(5). When aj is increased 
above this value, a very fast growth of the oscillation amplitude was observed. In particular, for aj = 0.333 the 
oscillation amplitude of the average z— projection of the element magnetization mf (the in-plane projection along 
the short ellipse axis, nearly perpendicular to the external field) was already A z rj 0.7 (see Fig. [IJ}). The regular 
oscillation regime exists in a extremely narrow a,/-region 0.331 < aj < 0.355. In this regime the spectrum of 
m| v -oscillations contains only a single very narrow peak which frequency is f± ~ 10.37 GHz for aj = 0.333 and 
decreases down to fx w 9.14 GHz for aj = 0.350. We point out that already for a,j — 0.333 the steady-state 
oscillation frequency is well below the frequency of the homogeneous FMR-mode /pw Rl although the distribution 
of the oscillation power is almost perfectly homogeneous (Fig. 2J), left panel). The reason is that even for this 
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aj-value - which is very close to the oscillation onset threshold a c J - the oscillation amplitude is so large that the 

formula /fmr~ (t/^ 71 ") • {Ho(Hq + AttMs)) 1 ^ 2 used to calculate /p 1 ^ (where small deviations of magnetization from 
its equilibrium orientation are assumed) is not applicable. 

As mentioned above, the regular oscillation regime with a well defined limiting cycles for the average system 
magnetization (Fig. right panel) persists up to Oj h w 0.355. Further small increment of the current strength leads 
to an abrupt transition to a quasichaotic magnetization motion (whether it is truly chaotic in the sense of non-linear 
dynamics should be investigated separately; see also our brief discussion in Ref. l26f) . An example of the magnetization 
trajectories in this quasichaotic regime is given in Fig. 

For currents not much larger than the transition value a^ h the oscillation spectra of the average magnetization 
components still exhibit a well defined structure with several peaks corresponding to several oscillation cigcnmodes 
with different spatial distribution of the oscillation power. An example of the such a spectrum with spatial maps 
attributed to different spectral peaks is given in Fig. ^ (left panel). In- plane spatial distributions of the oscillation 
power were calculated from the time-dependent trajectories of the discretization cell moments as described in Ref. 
1271 The detailed analysis of these spatial structures is beyound the scope of this paper. Here we would like only to 
mention that the asymmetry of the spatial distributions for the second and third modes is due to the presence of the 
Oersted field. Although this field is much smaller than the external field, its influence for this particular case is visible 
due to a relatively large amplitude of the magnetization oscillation. 

Further increment of the current strength (aj value) results in the gradual decline of the oscillation power for both 
in-plane magnetization components, accompanied by the decrease of the oscillation frequency as a function of aj. This 
behaviour is demonstrated in Fig. |3 where we present the maps of m!^f- and ra^ v -oscillation power as functions of 
the frequency (vertical axis) and aj (horizontal axis) together with the magnetoresistance oscillation power computed 
according to Eq. J2Jl. Characteristic for this dynamical regime is a multi-domain magnetization configuration whereby 
the average domain size decreases with increasing aj. Because the average magnetization is still directed along the 
external field, a formation of such a chaotic domain structure with smaller and smaller domains leads on average to 
the increase of the stray field in the opposite direction, which might be one of the reasons why the oscillation frequency 
decreases with aj. 

It is important to note, that in contrast to m| v , the oscillation power of the ra^ v -component (parallel to the long 
ellipse axis) is concentrated in two frequency regions. First, there is a contribution from fast m^ v oscillations at 
the frequency approximately twice as high as for m| v oscillations. The reason for this relation between m x and m z 
oscillation frequencies is that for the motion type under study - the magnetization is precessing around (approximately) 
the Ox-direction - the m^-componcnt moves back and forth twice during one complete oscillation cycle (see, e.g., right 
graphs in Fig. ^3 and c). Second, there exist a substantial low- frequency contribution to the oscillation power (which 
growth with increasing a,) due to slow variations of m^ v (t) which are typical for a quasichaotic mo tion. A nalogous 
behaviour of the longitudinal magnetization component in similar situations was observed in Ref. I2fil28l2°l3f)l (see 
also the discussion below). 

Increasing aj above the value aj ~ 1.4 leads to the disappearance of the chaotic regime and to the establishing 
of a so-called 'out-of-plane' precession mode. The magnetization trajectory in this mode is an ellipse with the plane 
slightly tilted with respect to the nanoclemcnt plane. The oscillation frequency abruptly increases and is now for 
obvious reasons the same for m^ v and m| v -components (see regions marked by circles in Fig. [2Jl . The magnetization is 
precessing nearly coherent in this mode, what can be seen from a large absolute value of the out-of-plane magnetization 
projection: (m™ « 0.8 — 0.9). The absolute value of the oscillation power in this regime is at least an order of 
magnitude smaller than at the beginning of the quasichaotic motion due a relatively small precession amplitude. 

Comparison of these results with experimental data and numerical simulations of other groups will be made in 
Section HV| 



B. Magnetization dynamics for T = 0: influence of the exchange stiffness value 

Taking into account a relatively broad region of the exchange stiffness constants reported for bulk Co and thin 
Co films (see discussion at the beginning of this Section), we have studied how the decrease of A would affect the 
magnetization dynamics. From general considerations, we expect the magnetization configuration to become 'softer' 
for smaller exchange stiffness, so that the influence of a strongly inhomogeneous self-demagnetizing field of a flat 
elliptical nanoelement should lead to an equilibrium magnetization configuration which is less collinear compared to a 
system with larger A. This, in turn, should result in a broader distribution of the resonance fields for different regions 
of a nanoelement (different discretization cells in our micromagnetic model) , thus leading to a slower increase of the 
magnetization oscillation amplitude with increasing aj (above the oscillation threshold a"). In addition, for smaller 
exchange the transition from a regular to a quasichaotic oscillation regime is expected to occur earlier (for smaller 
aj)- 
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All these features were indeed observed in our simulations, whereby it has turned out that from the quantitative 
point of view the nanoelement with the geometry under study is quite sensitive already to moderate decrease of A. 
In Fig. 01 we display an example of the transition from a regular to a quasichaotic behaviour for the nanoelement 
with parameters identical to those given in the preceding subsection, but with A = 2.0 • 1CP 6 erg/cm. The oscillation 
threshold a c J — 0.3290(5) is almost the same as for the system with the reference parameter set. However, the 
spectral power of magnetization oscillations in a regular regime is substantially smaller than for a 'reference' system. 
In particular, when a,j exceeds the oscillation threshold value by Aaj = 0.0015 (aj = 0.331, see corresponding 
spectrum in Fig. and trajectories in Fig. |3p), the oscillation amplitude is significantly smaller than for the same 
situation for a 'reference' system (Fig. QJd). For this reason the oscillation frequency here (/ » 12.45 GHz) is much 

closer to /p^pthan for the peak shown in Fig. \I]p. 

The transition to a quasichaotic motion occurs, due to the reason explained above, for a somewhat smaller current 
value a} h = 0.338(1). Interestingly, the non-trivial mode structure appears for this system already before the transition 
(see spectra for aj = 0.334 and aj = 0.336 in Fig. [^Ji). In addition, for the system with smaller A only the first 
and third mode from of the eigenmode set identified for A = 3.0 • 10 -6 erg/cm are visible in the power spectrum of 
magnetization oscillations just above the transition to chaos. The analysis of trajectories for individual discretization 
cells reveals that the second mode (with the intermediate frequency) is still present. However, the spatial structure of 
this mode for this smaller exchange value is nearly symmetric (compare oscillation power maps in Fig. ^ and|3t), and 
the mode does not manifest itself in the spectrum of the average magnetization dynamics, because the magnetization 
in the left and right mode localization regions oscillates in opposite phases. 

Analysing the influence of exchange stiffness value on the overall spectral picture in the chaotic regime, we can 
see that the decrease of A increases the upper boundary of the current region where the quasichaotic regime exists 
(compare Fig. [SJi andlSJa). In addition, it can be seen that the decrease of A down to A = 2.0 ■ 10 -6 erg/cm leads 
to a nearly complete disappearance of the 'out-of-plane' precession region. This means that for large spin torque 
values this exchange constant is not high enough to ensure the collinear magnetization structure necessary for the 
very existence of the 'out-of-plane' precession mode. 

C. Magnetization dynamics for T = 0: influence of the random crystal grain anisotropy 

It is well known that thin sputtered Co films have a polycrystalline structure which in the simplest case can be 
treated as absolutely random and characterised by a single parameter - the average grain size (D) . In all the simulations 
which results are shown below we have used (D) = 10 nm. To study a nanoelement with such a random polycrystalline 
structure we have applied a standard procedure included into the MicroMagus^i package. The procedure starts with 
the random placement of 'growth centers' within the simulated nanoelement, whereby the number of growth centers 
is determined from the nanoelement size and the average grain size (D). Afterwards the crystal grains are 'build' 
simulating a simple isotropic 2D growth which starts simultaneously at all growth centers and is terminated at the 
locations where the adjacent grains touch each other. The routine stops as soon as the whole element area is covered 
by the grains. Finally a random direction of the anisotropy axes is assigned to each grain, which means that (i) for 
all discretization cells within a given grain the anisotropy axes directions are the same, but (ii) the corresponding 
directions for two different grains are not correlated. 

Thin Co films may have a polycrystalline structure with grains having either fee or hep crystal lattice type. Co 
films with thicknesses about several nanometers epitaxially grown on Cu substrates possess the fee crystallographic 
structure (see, e.g., review 3 ^ and references therein), so it is likely that Co films sputtered on the same substrate (as 
those used in Ref. H4T) also have predominantly the fee structure (see also results and discussion in Ref. l3Sh . Hence 
we begin this subsection with the analysis of simulation results for the corresponding case. 

1. Magnetization dynamics for a nanoelement with fee crystal grains 

At room temperatures the fee crystallographic modification of Co is unstable for bulk crystals, so that measurements 
of the anisotropy constant can be performed only on thin films which fee structure is stabilized by, e.g., a suitable 
substrate choice (like Cu with the Ill-orientation of the substrate surface). This circumstance leads to considerable 
difficulties by the anisotropy constant determination, which could be performed mainly on epitaxially grown films (see, 
e.g., Ref. l39l) . The values of the first cubic anisotropy constant obtained for thin Co films grown at different conditions 
on different substrates are within the range Kf nh = (5 — 8) • 10 5 erg/cm 3 , so we have used the value XJ ub = 6.0 ■ 10 5 
erg/cm 3 in our simulations. Varying Kf uh within the region cited above did not lead to any significant changes of the 
results. 
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From this anisotropy value and the standard saturation magnetization of Co M$ ~ 1400 G (which is almost the 
same for fee and hep phases) one can determine the reduced anisotropy constant — 2i"T 1 /Mj, also sometimes referred 
to as a 'quality factor' Q. The corresponding value (3 ~ 0.6 is not small (and is even larger if we substitute into the 
expression for j3 the reduced Ms value reported in Ref. Il4|) . which means that this random anisotropy may cause 
substantial inhomogeneities of magnetization structure in a polycrystalline sample. In nanocrystalline thin films, 
however, where the grain size is of the same order of magnitude as the exchange length, the influence of the random 
anisotropy is strongly diminished, because its effect is 'averaged out' due to uncorrelated anisotropy directions in 
adjacent crystallites^. 

Nevertheless in the nanoelement under study (with the average grain size {D) = 10 nm) the random cubic anisotropy 
leads to a noticeable deviations of the magnetization structure from the configuration obtained for an ideal nanoelement 
(without the polycrystalline structure) and thus results in an observable spread of resonance fields for different 
crystallites. Qualitatively this leads to the same effect as the decrease of the exchange constant: after passing the 
threshold value a c j = 0.337(2) the amplitude of the SPC-driven precession of the average magnetization increases 
much slower, than in the absence of the random anisotropy. In particular, for aj = 0.340 we can still observe a really 
small-angle homogeneous precession with the frequency /hom ~ 13.6 GHz which is virtually equal to /fmr(~ 1^.7 
GHz) - see Fig. Hd. 

In contrast to the decrease of the exchange constant, which leads to systematic deviations of the magnetization 
configuration from a homogeneously magnetized state under the influence of the self-demagnetizing field, a crystal 
grain anisotropy leads to random deviations of the magnetization. For this reason the mode structure after the 
transition to the quasichaotic precession (Fig. 0J;) looks entirely different when compared with the case of the reduced 
exchange constant (Fig. |3J:). And finally, the analysis of the oscillation spectra in the quasichaotic region of currents 
for the nanoelement with a random polycrystalline structure (Fig. Et) shows that the presence of such an anisotropy 
slightly expands the existence region of quasichaotic oscillations (compare with Fig. [SJi) and significantly decreases 
the spectral power of 'out-of-plane' oscillations. The latter effect is due to the same reason as for the case of a reduced 
exchange - random anisotropy disturbs the homogeneous magnetization configuration required for the existence of an 
'out-of-plane' precession. 

2. Magnetization dynamics for a nanoelement with hep crystal grains 

It has been shown experimentally^ that although crystallites of thin Co films sputtered on Cu/111 substrates have 
mostly the fee structure, there exist a non-negligible fraction of grains possessing a 'normal' hep structure of bulk 
Co crystals. For this reason we have also studied the magnetization dynamics of a polycrystalline nanoelement with 
hep-grams having the same average size (D) = 10 nm. 

Magnetic anisotropy of the hep modification of Co can be treated for our purposes as a uniaxial one with only the 
first anisotropy constant Kf ub = 4.5 • 10 6 erg/cm 3 (see, e.g., Ref. Ell) taken into account. The reduced anisotropy 
constant (3 w 4.6 calculated from this K\ value and the same saturation magnetization Ms w 1400 G is so large 
that the random anisotropy plays a crucial role even in a sample with such small crystallites. This can be seen 
very well in Fig. where we show maps of the oscillation power for polycrystalline nanoelements with hep grains. 
The four maps (a)-(d) were computed for 'samples' with identical macroscopic parameters ('reference' parameter set 
and one and the same value of the random uniaxial anisotropy of crystallites Kf uh = 4.5 • 10 6 erg/cm 3 ), but using 
different realizations of a random polycrystalline structure. The most striking feature of the spectra shown in Fig. 
is that they are qualitatively different for different 'samples'. In particular, we can recognize both (i) a large spread 
of the regular oscillation frequencies (narrow peaks for small aj lie in a wide range « 10 — 20 GHz) and (ii) very 
different width and aj-dependencies of broad spectral bands in the region of a quasichaotic behaviour. In addition, 
we observe large frequency jumps by the transition from a regular to a quasichaotic regime (see, e.g., realizations (a) 
and (b)). This is in strong contrast to the behaviour of a nanoelement with an absent or moderate anisotropy, where 
the regular oscillation frequency before the transition to a quasichaotic regime was 'inherited' by a main spectral peak 
of quasichaotic oscillations just after the transition. 

Further insight into the magnetization dynamics of a polycrystalline nanoelement with hep grains can be gained 
by the analysis of the magnetization trajectories and spatial distribution of the oscillation power. 3D trajectories of 
the average system magnetization for the 'sample' marked as (b) in Fig. El at the aj values before (aj = 0.50) and 
immediately after (aj = 0.54) the transition from a regular to a quasichaotic regime are shown in Fig. [7| together 
with m^ v and mj v power spectra and spatial maps of m^ v and m| v oscillation power. 

First of all we note that due to a strong influence of the hep anisotropy with randomly oriented anisotropy axes 
in various grains, the average nanoelement magnetization does not precess anymore approximately around the Ox 
axis (long ellipse axis) even in the regular regime. This is clearly demonstrated in Fig. EH: the limiting cycle in 
the form of a narrow bent torus is strongly tilted with respect to the Ox-axis. Such an asymmetric trajectory lead 
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to the appearance of a strong spectral peak at the basic precession frequency also in the oscillation spectrum of the 
longitudinal magnetization component m^ v (compare with spectral maps for an element without a random anisotropy 
in Fig. Ek,b). 

The next important feature is a strongly irregular spatial distribution of the oscillation power which is especially 
pronounced for the m^ v -component (grey-scale map in Fig. 03). It is evident that this power distribution is determined 
by the random polycrystalline structure; this circumstance is again in a strong contrast to the case of an absent or 
moderate fee anisotropy. By the transition to a quasichaotic regime this spatial distribution changes drastically 
(compare corresponding maps in Fig. 03 and 03) thus causing a relatively large frequency jump by this transition. 
For example, the difference between the positions of the peaks corresponding to basic oscillation frequencies in Fig. 
03 and 05 is » 3.5 GHz. 

Further discussion of these results is again postponed till the next Section. 

D. Magnetization dynamics for finite temperature: influence of thermal fluctuations 

Simulations of the magnetization dynamics taking into account thermal fluctuations were performed by adding the 
'fluctuation field' H fl to the total effective field in the LLG equation of motion. We have used components of H fl with 
zero means and (5-correlated in space and time 

(HlMH^it)) =2D5^S id 8(t) (3) 

where the noise power D — (A/1 + A 2 ) • (fcT/7//) depends on the system temperature T, damping A and the total 
magnetic moment of a discretization cell [i = MgAV. In principle it is known that the non-trivial correlation 
properties of the random noise appear in numeric al micr omagnetics due to the discretization of a formally continuous 
problem. We refer the interested reader to Ref. I27l42l for the detailed discussion about these properties and their 
possible influence on the magnetization dynamics. Here we only note that we have checked the issues mentioned in 
Ref. I27l42l to ensure that the approximation J3) is sufficiently good for our situation. 

Another problem arising by numerical simulations at finite temperatures is whether one should change the macro- 
scopic magnetic parameters of the material under study when doing such simulations (see Ref. |43j and citations 
therein for the recent discussion). Taking into account that we use fairly small discretization cells and simulate the 
system behaviour at room temperature (T — 300 K) which is well below the Curie temperature of common magnetic 
materials, we have left the macroscopic parameters unchanged. However, we note that the arguments presented in 
Ref. ^3 make further studies in this direction worthwhile. 

We also point out that the inclusion of the fluctuation field H fl with correlation properties given by © converts 
the ordinary differential LLG-equation into a stochastic one. Due to random functions appearing on the right-hand 
side of such equations numerical methods for their solution have, in general, a much lower accuracy order than their 
counterparts for ordinary differential equations 44 . For this reason a much smaller time step is required to integrate 
the LLG-equation with the prescribed accuracy at T > 0. In our particular case, for T — 300 K and discretization 
cells with a relatively small size (note that the noise power in Eq. J3J is inversely proportional to the cell volume), the 
time step had to be decreased by more than order of magnitude. Hence for T — 300 K we have performed simulations 
only during the physical time of « 50 ns (after establishing the steady state oscillation regime). This circumstance 
explains why the statistical errors in spectra shown in Fig. are considerably larger than for those analysed in 
previous subsections (for T = 0), especially for a quasichaotic oscillation regime. 

Example results of our simulations including thermal fluctuations are displayed in Fig. where we present the 
analysis of the magnetization dynamics of a polycrystalline nanoelement with a 'reference' parameter set and fee Co 
grains. 

Comparing these data with the results shown in Fig. and 0; we can see, first of all, that thermal fluctuations 
have a qualitative influence on the regular precession regime (the region of regular oscillations is marked also in Fig. 
|SJ. On the one hand, the aj-region where corresponding oscillations exist, is much broader than for T = 0: already 
for aj = 0.2 the m^ v oscillation power is twice as large as in the absence of a spin current. On the other hand, the 
width of corresponding spectral lines drastically increases. For T — the linewidth of m| v -oscillation spectra for a 
regular precession was always Af < 20 MHz, whereby this number represents an upper limit posed on Af by the finite 
duration of our simulations. Inclusion of a thermal noise leads to strong aj-dependent broadening of the spectral lines: 
at room temperature in the regular regime Af decreases with aj starting from the value Af 1.5(±0.2) GHz for the 
spectrum of spontaneous oscillations (without a spin current, aj = 0) down to Af w 0.4(±0.1) GHz for aj = 0.3. 

The second effect of the thermal noise is the strong decrease of the 'out-of-plane' oscillation power. Comparison 
of Fig. 03 and 0: shows that in the presence of thermal fluctuations this power is smaller by more than one order 
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of magnitude. The reason is the same as explained above by the comparison of results for a 'minimal model' and its 
modifications presented in Fig. [5] 

And finally we note that according to our simulations, thermal noise does not change the frequencies of spectral 
peaks in a region where a quasichaotic precession regime is well established (starting from aj « 0.38 for the system 
studied here). Comparison of corresponding oscillation spectra reveals that the frequencies of both peaks remains the 
same in frames of statistical errors. 



IV. DISCUSSION 
A. Comparison with numerical simulation results of other groups 

Full-scale micromagnetic simulations of a magnetization dynamics driven by a spin-polarized current have been 
performed, up to our knowledge, only in a few papers, which are partially mentioned in the Introduction. 

The first systematic study of the SPC-driven magnetization dynamics was performed in Ref. where a noise 
induced in a double-layer system with the geometry corresponding to a CPP spin valve head was investigated. Several 
features of the spin-current induced precession observed in Ref. |28| (fast increase of the magnetization oscillation power 
after the appearance of a steady state precession, linear dependence of the current threshold value on the damping 
constant A, uniform magnetization precession for low current values and chaotic magnetization motion for large 
currents, significant low-frequency contribution to the oscillation power for large currents) are qualitatively similar to 
our results. However, a quantitative comparison is not possible, because geometry and material parameters used in 
Ref. |2^ are very different from ours. 

The steady-state magnetization precession in an elliptical single-layer nanoelement 'made of Permalloy (NisiFeig) 
was studied for the first time in the short report^. A polycrystalline structure of a nanoelement was neglected taking 
into account an extremely small anisotropy of Py grains. Both regular and quasichaotic precession regimes were 
found in Ref. for an element with lateral sizes 100 x 50 nm. The increase of the oscillation amplitude beyound the 
threshold of the steady-state precession onset observed in Ref. [2!| is somewhat slower than in our simulations. The 
current region (normalized on the precession onset value) where a regular precession regime exists in Ref. is broader 
than for our nanoelements without a polycrystalline structure. Both quantitative discrepancies can be explained by 
the smaller exchange constant A = 1 x 10~ 6 of Permalloy (compared to our system) and a smaller element size used 
in Ref. |U 

Simulations most closely related to our studies were done in Ref. I30L where the magnetization dynamics of a 
single-layer elliptical element with the same geometry as used by us was performed. Lee et al»22i did not study how 
the magnetisation dynamics depends on the nanoelement magnetic parameters and did not carry out a systematic 
analysis of the oscillation spectra, but it is nevertheless possible to compare some results from Ref. |33 with our 
findings. 

First of all we note that Lee et al. simulated a nanoelement with the standard saturation magnetization value of 
the bulk Co Ms = 1400 G. Due to this Ms-value they could obtain a satisfactory agreement between the oscillation 
frequency in a regular precession regime (called a " small-amplitude signal" in Ref. |3(3 and fH) obtained in their 
simulations and experimentally^. For example, for the external field Hq = 2 kOe simulations performed in Ref. l30l 
produce / s j m rj 16.5 GHz (see Fig. 3b in Ref. 30), whereby experiment of the Cornell group gives— f cxp « 16 GHz. 
This agreement is a natural consequence of a higher Ms value used in in Ref. |3(] compared to our simulations and the 
Ms-value measured for a sample used in Ref. [Tip the Kittel FMR-frequency obtained for a thin film with Ms = 1400 
G at Hq — 2000 Oe is /pjIr ~ VI .h GHz. An slightly smaller value simulated in Ref. 30 is most probably due a 
slightly non- linear precession character (see our discussion of results presented in our Fig. 0. 

We have also performed simulations with Ms = 1400 G (results are not shown here) and have obtained for a 
nanoelement with the same geometry a regular precession with the frequency / S i m « 17.2 GHz, which is very close 
both to the value obtained in Ref. y(J and that calculated from the Kittel formula. However, the quasichaotic regime 
for an element with such a high saturation magnetization starts by our simulations at the frequency « 16 GHz what 
is approximately 3 times larger than the experimentally observed value (see Fig. lc and If in Ref. IT4T) . In this 
connection we note, that although the frequencies for a quasichaotic regime are not given in Ref. l30t the transition 
scenario from a reg ular to a quasichaotic precession shown in Fig. 3c from Ref. M to similar to ours (see, e.g., Fig. 
n in this paper) , so that we expect that the frequencies at the beginning of a quasichaotic motion for simulations 
performed in Ref. |3(3 are also approximately as high as those for regular oscillations. 

In addition, the relation of the maximal to the minimal current of the range where the quasichaotic regime exists for 
a nanoelement with Ms = 1400 G is also much larger than the value / max /^min ~ 3 deduced from the data presented 
in Fig. If of Ref. [14]. The relation J ma x/-^min increases, because the orientation (on average) of the magnetization 
along the external field is preserved up to higher current values for a nanoelement with larger Ms- Summarizing, the 
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usage of a larger saturation magnetization changes the overall agreement between the experiment and simulation for 
the worse. 

It is also instructive to compare our results with the simulations performed in the macrospin approximation in Ref. 
Phi itself. First of all we note that by performing such simulations one has to introduce the artificial anisotropy field 
H an (directed in this case along the long axis of the ellipsoidal element) in order to reproduce already the quasistatic 
hysteresis loop. The magnitude of this field (or the corresponding anisotropy constant) is usually chosen to obtain 
the coercivity similar to those measured experimentally. From the physical point of view this field is intended to 
describe the average magnetocrystalline anisotropy of a polycrystalline element and/or its shape anisotropy. The first 
contribution is negligibly small in samples without a specially induced in-plane grain texture (what is probably absent 
for sputtered Co films). The second contribution - the shape anisotropy of a thin elliptical nanoelement - can not 
be described adequately by a uniaxial anisotropy term alone, but is in fact a superposition of a strong 'easy plane' 
anisotropy due to the demagnetizing field of a thin film and an additional in-plane anisotropy arising due to a non- 
circular lateral shape of a nanoelement. This is one of the reasons why one can expect only a qualitative information 
from such a macrospin approximation, even if no domain formation is expected during the magnetization reversal. 

From this qualitative point of view macrospin simulations presented in Ref. demonstrate two interesting features 
(see Fig. 3 in Ref. 14): (i) a frequency jump by the transition from a small-angle to a large-angle precession and (ii) 
an out-of-plane precession mode with the frequency increasing with the current strength. 

The first feature is temptingly similar to the frequency jump found for the sample 1 from Ref. ^3 when the current 
was increased (see Fig. Id and If) and hence was proposed by Kiselev et al. as a possible explanation for this exper- 
imental observation. Our micromagnctic simulations reveal, however, that in a full micromagnetic model this jump 
does not occur by the transition from a regular to a quasichaotic precession regime. The reason for this discrepancy 
between the macrospin and full micromagnetic model is the following. In the macrospin approximation the abrupt 
transition between the small-angle (ellipsoidal) and a large angle ('butterfly' or 'shell') precession is accompanied 
by the strong decrease of the oscillation frequency because the length of the 'butterfly' trajectory after the transi- 
tion is much larger than the length of the elliptical trajectory before it. Full-scale micromagnetic simulations reveal 
that the regular precession regime exists really up to nearly maximal amplitudes of the magnetization oscillations. 
This large amplitudes can not be significantly increased by the quasichaotic regime and thus the precession in this 
regime 'inherits' the frequency of regular oscillations, and the frequency jump by this transition is absent. Hence the 
abovementioned transition can not be the reason for the frequency jump observed in a real system. 

The second feature obtained in the macrospin approximation - the 'out-of-plane' precession regime - is not observed 
experimentally 14 . Indeed, our full-scale simulations show that this precession mode is relatively unstable due to a 
slightly non-collinear magnetization configuration already for the 'minimal' model, where neither the polycrystalline 
structure nor thermal fluctuations are taken into account. Either of these two factors disturbs further the homogeneous 
magnetization structure, thus nearly wiping out the 'out-of-plane' mode, which intensity is then probably below the 
observable threshold. 



B. Comparison with experimental data 

As it was mentioned in the Introduction, the quality of experimental data presented in Ref. ^] is sufficiently high 
to allow a quantitative theoretical analysis, but one should of course begin with a qualitative comparison between 
experimental and simulated results. 

First we note that the presence of several spectral bands shown in Fig. If in Ref. ^] is due to the superposition 
of signals arising from the oscillations of the perpendicular and longitudinal magnetization components. This can 
be clearly seen from Eq. J2J an d Fig. where the magnetoresistance oscillation spectrum obtained using @ 
is presented. Oscillations of the longitudinal (m| v ) magnetization component produce the signal with the basic 
(precession) frequency, whereby the magnetization component perpendicular to the long axis of the ellipsoid (rn^) is 
responsible for the second harmonic generation: As explained above by the analysis of the 'minimal model' spectra 
( section IIII A|) . for this type of motion the m^ v -component oscillates with the frequency roughly twice the precession 
frequency. A weak third harmonics is present in the spectrum of m^ v -oscillations due to a non-linear precession 
character for this regime. A quantitative comparison of the basic frequency and harmonic intensities is not possible, 
because the relation between corresponding intensities strongly depends on the angle between the external field and 
the ellipsoid axis, which is not known exactly. 

The interpretation outlined above is strongly supported by the experimental observation that the basic frequency 
band disappeared when the external field was aligned as exact as possible along the long axis of the elliptical element^, 
as it should be according to Eq. J2J with e h z = 0. 

The next intriguing feature of the magnetization dynamics measured in Ref. ^3 is t he existence of two oscillation 
regimes - for low (for 2.0 — 2.4 mA) and high (from « 2.6 mA up to 6.5 mA) current values. The transition between 
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these two regimes is accompanied by the drastic increase of the microwave power and by the frequency jump from 
rj 16 GHz in a low-amplitude regime down to « 7 GHz for large-amplitude oscillations. These features are quite 
similar to those obtained for the macrospin model (see corresponding simulations in Ref. IT3) . However, as explained 
above, they disappear when this transition is studied using a full micromagnetic model with moderate or absent 
random grain anisotropy of a polycrystallinc elliptical clement: First, the frequency during this transition remains 
nearly unchanged, and second, the magnitude of the simulated microwave power before the transition is at least as 
high (and for T = even higher) as immediately after it (this latter feature is also due to a fact that the amplitude 
in a regular regime just before the transition is by no means small, as it is the case in the macrospin model). 

For this reasons another explanation of the frequency and amplitude relation in these two regimes is required. Such 
an explanation can be given assuming that the nanoelement studied in Ref. ^| had at least a small (but non- negligible) 
fraction of hep Co grains (evidences about the coexistence of fee and hep phases in this Co films can be found, e.g., 
in Ref. l38|) . As shown in Fig. polycrystalline 'samples' consisting of hep Co demonstrate a large variety of possible 
spectral maps which strongly depend on the realization of a polycrystallinc sample structure. As explained in section 
IIII C 21 this is a natural consequence of a relatively strong anisotropy of the hep grains, which lead to very different 
magnetization dynamics for garins with different directions of the anisotropy axes. This, in turn, may lead to a spatial 
'migration' of the oscillation power by increasing the current strength. All these circumstances can result in a single 
(Fig. [U)) or multiply (Fig. ®l) frequency jumps when increasing current and quite different dependencies of the 
spectral maxima frequency on the current strength (compare the behaviour of broad spectral bands in, e.g., Fig. EJ; 
and d). The magnitude of the frequency jump can be quite large: we have observed frequency changes as large as 10 
GHz (in principle even larger jumps are possible: the maximal oscillation frequency for a hep Co grain with Ms of a 
bulk Co and the anisotropy axis aligned along the external field is / max » 38 GHz) . In this connection it is interesting 
to note that the regular precession frequency / rog w 17 GHz obtained for the sample which spectral map is shown in 
Fig. EJd nearly coincides with that measured in the 'small-amplitude' regime in Ref. Il4l 

The interpretation proposed above could also explain why the magnitude of the microwave power drastically in- 
creases after the transition to a quasichaotic behaviour: it might happen that only a small fraction of a sample (e.g., 
only hep grains with a favourable orientation of the anisotropy axes) oscillates in a regular ('small-amplitude') regime, 
whereby in a quasichaotic regime one observes the oscillations of a major part of a nanoelement occupied by the fee 
nanocrystals. Our hypothesis is supported further by the observation that the 'small-amplitude' regime was found 
not for all experimentally studied samples^. 

We would like also to add a few comments concerning the 'out-of-plane' precession mode which was not observed 
experimentally, despite being predicted by a macrospin model. As discussed in subsections IIII C II IIII Dl and IIV Al 
above, the absence of this mode in the measured spectra can be explained by the joint influence of a polycrystalline 
sample structure and thermal fluctuations: both these factors destroy the homogeneous magnetization structure thus 
decreasing the amplitude of the 'out-of-plane' mode. In this connection it is worth noting that for some /icp-'samples' 
we have observed a relatively weak spectral band corresponding to 'out-of-plane' oscillations (this mode could be 
identified due to its very characteristic frequency dependence on the current strength similar to that shown in Fig. 
3b in Ref. IbH) . We have found out that this mode was localized in a small region of a nanoelement, where the 
orientations of the grain anisotropy axes were nearly parallel, thus favouring the collinear magnetization structure 
during the precession. The importance of the polycrystalline structure for the suppression of the 'out-of-plane' 
precession could be tested on samples made of materials with low crystal grain anisotropy (like Permalloy). For this 
samples the 'out-of-plane' mode should be observable, because the nearly absent magnetic anisotropy of Py could 
not destroy the homogeneous magnetization configuration required for the existence of the 'out-of-plane' mode. This 
conclusion seems to be supported by the recent experimental observations^. 

As it can be seen from the previous discussion, the extended micromagnetic model which takes into account a 
polycrystalline structure of a nanoelement and thermal fluctuations can explain qualitatively all most important 
features of the experimental data reported up to nowi^. However, from the quantitative point of view, several 
discrepancies do exist. 

The most important one is the difference between the simulated and measured frequencies in the quasichaotic 
regime: at the beginning of this regime the measured frequency is / mos ~ 6.5 GHz (see Fig. lc in Ref. \l4) , whereas 
our simulations for elements without a polycrystalline structure and for a polycrystalline fcc-Co give / s j m w 9 — 10 
GHz for a j values where the transition to a quasichaotic regime is completed (see, e.g., Fig. {lfc and spectrum for 
a j = 0.38 in Fig. [SJ. We consider this discrepancy as the most serious one because the frequency is the most reliable 
spectral characteristic obtained both in measurements and simulations: whereas other parameters like the line width 
and intensity can be artificially altered by many factors, the peak frequency is an inherent signal feature which is 
usually reproduced with an accuracy allowing a direct comparison between simulations and experiment. 

By the discussion of the influence of various factors which could account for the frequency difference pointed out 
above we would like to begin with the reasons having no direct relation to the spin injection, i.e., with the factors 
which are not taken into account in our 'standard' micromagnetic model. 
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(i) Influence of the underlayer. In our treatment we have considered a single-layer elliptical nanoelement, which 
means that the magnetodipolar interaction with the Co underlayer present in a real experiment^* was not included. In 
principle, such an interaction could indeed decrease the oscillation frequency, especially taking into account that the 
underlayer could be also partially structured by the ion milling used in Ref. 14 to produce an elliptical nanoelement 
from the upper Co layer. The elliptical pillar sticking out from the underlayer after such a milling would be magnetized 
in the external field direction thus producing the stray field in the opposite direction. This latter field would decrease 
(on average) the total field acting on the elliptical nanoelement under study thus decreasing its precession frequency. 

We are not able to simulate this effect quantitatively because the milling depth for the underlayer is not known 
exactly. However, by creating the columnar structure studied in Ref. it was aimed to disturb the lower Co layer 
as little as possible, so that the milling depth through the underlayer should be maximal several nanometers*^. The 
stray field from a pillar with such a small height acting on the element placed 10 nm above it would be concentrated 
mainly in small regions near those edges of the nanoelement which are perpendicular to the external field (or the long 
ellipsoid axis). Thus such a field could presumably not lead to a significant decrease of the precession frequency in 
the quasichaotic regime, where oscillations of the inner nanoelement regions play a major role. 

To check this conclusion, we have performed simulations of a system consisting of an extended underlayer with a 
thickness 10 nm, an elliptical pillar with the thickness 5 nm and lateral sizes 130 x 70 nm placed on the top of it 
and a nanoelement with the same lateral sizes located 10 nm (thickness of the Cu layer between the two Co layers 
in Ref. fl4|) above the pillar. Periodic in-plane boundary conditions were applied to eliminate the effect of the stray 
field from a relatively thick extended underlayer. We have found (detailed results will be presented elsewhere) that 
the precession frequency both in the regular and in the quasichaotic regime was not noticeably influenced by the 
underlayer with a pillar, as expected according to the considerations above. 

(ii) Grain texture for an element with a hep polycrystalline anisotropy. The already studied influence of the poly- 
crystalline structure for a Co nanoelement with hep grains could be greatly enhanced if the element would possess a 
significant grain texture. In particular, if there would exist an element region with the lateral size much larger than 
the element thickness where the grain anisotropy axes are approximately parallel to each other and perpendicular 
to the element plane, then the oscillation frequency for a homogeneous precession mode of such a region would be 
f^ R = (7/271-) • [Ho (H + (4?r - (3)M S )} 1/2 (for (3 < 4?r), where (3 = 2K/M% is the reduced anisotropy constant intro- 
duced above. For a hep Co with Ms = 1400 G and K — 4.5 x 10 6 erg/cm 3 this means the reduction of the precession 

frequency by a factor /fmr//fmr ~ 1-2- If we assume that the quasichaotic precession frequency / c h scales with 
Ms and (3 at least approximately in the same way as the FMR-frequency, the influence of the /icp-anisotropy could 
decrease / c h from 10 GHz down to « 8.2 GHz which would significantly improve the agreement with the experiment. 
For a Co film with the reduced Ms ~ 800 G the anisotropy (3 ~ 14 would be even somewhat larger than the shape 
anisotropy so that the precession with the frequency close to the free spin value fo = (7/2-71-) • Ho ~ 5.5 GHz and 

lower (taking into account that / c h < /fmr) cou ld be expected. Unfortunately, as stated above, a highly textured 
polycrystalline structure should exist in the sample under study in order to make these arguments applicable. This 
is unlikely for a sputtered Co film (we note in addition, that such an area would probably lead to an appearance 
of a significant peak at the corresponding frequency in a regular regime). However, this line of the argumentation 
clearly shows that a detailed characterization of a polycrystalline sample structure is necessary when one performs 
measurements on materials with moderate to high magnetocrystalline anisotropy. 

(iii) Surface roughness. In a thin film element with a rough surface a demagnetizing field due to this roughness 
would exist, what can also decrease a precession frequency. Quantitative studies of the roughness effect are possible 
only when the corresponding information about the roughness parameters is available. Here we can only point out 
that a substantial roughness with a characteristic wavelength of the same order as the element thickness is necessary 
to create an appreciable demagnetizing field inside a nanoelement. 

Next we consider in our discussion the possible factors connected to the magnetization dynamics of a system driven 
by a spin-polarized current (SPC). 

(iv) Effective field induced by the spin-injection. The existence of such a field was extensively discussed theoretically 
(see references in the Introduction), but only recently this field could be measured in a direct experiment^ using a 
small asymmetry of the hysteresis loop in a multilayer system in the presence of a SPC. For a system very similar 
to that studied in Ref. [lj it was shown that the magnitude of a SPC-induced effective field can be calculated as 
-ffspc = bj ■ j, where j is the current density and the coefficient bj was found to be bj « 1.5 x 10~ 7 Oe x cm 2 j A. For 
a typical current value 7 = 3 mA and lateral sizes of the nanoelement under study 130 x 70 nm the current density 
is j » 4 • 10 7 A/cm 2 . Hence a SPC-induced field would be -ffspc ~ 6 Oe which is negligibly small. 

(v) Non-trivial angular dependence of the Slonczewski torque magnitude aj and configuration- dependent damping 
A. As it was discussed in Sec. [H] these factors should be considered on equal footing, because the spin-induced torque 
and precession damping in multilayer systems subject to a spin-polarized current are controlled by the processes of the 
same physical nature. Here we would like only to mention, that from a qualitative point of view the enhancement of 
the precession damping for large SPC strengths and/or for large angles between the magnetization of adjacent layers 
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(see Ref. in Sec. |HJ could indeed decrease the precession frequency in a quasichaotic regime which is characterized 
by strong inhomogeneities of the magnetization structure. The quantitative study of the corresponding effect is the 
subject of a future work. 



V. CONCLUSION 



In conclusion, we have performed a systematic study of the magnetization dynamics of a thin single-layer elliptical 
nanoelement driven by a spin-polarized current flowing perpendicular to its plane. Analyzing the dependence of this 
dynamics on various magnetic material parameters, we have demonstrated that an extended micromagnetic model 
which takes into account a polycrystalline structure of this element and the influence of thermal fluctuations can in 
principle explain qualitatively almost all features observed experimentally. At the same time we have shown that 
an apparent qualitative agreement between the macrospin model and some experimental findings is due to artificial 
features of the macrospin dynamics, which are absent in the full-scale micromagnetic model. An important issue 
concerning the quantitative disagreement between the precession frequencies measured experimentally and simulated 
numerically for the quasichaotic regime still requires further clarification. Analyzing various possible causes of this 
discrepancy we have shown that the most likely candidate to resolve this difficulty is the refined model of the spin- 
injection induced dynamics which would include the non-trivial angular dependence of the spin-torque term and the 
dependence of the precession damping on the spin-polarized current strength and the magnetization configuration of 
a multilayer system. 
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FIG. 1: Magnetization dynamics for an elliptical element with a 'reference' parameter set near the oscillation onset threshold. 
Panel (a): power spectra of the mjf-component for various aj values as indicated in the figure. Panel (b): magnetization 
dynamics for aj = 0.333 represented by the enlarged mjf-oscillation spectrum and the projection of the average magnetization 
trajectory on the coordinate plane Oxz (coinciding with the nanoelement plane). Panel (c): the same as on (b) for aj — 0.360 
- immediately after the transition to a quasichaotic oscillation regime. Grey-scale maps represent spatial distributions of the 
mf oscillation power for corresponding spectral peaks. 
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FIG. 2: (Color on-line) Spectral power of magnetization oscillations for an elliptical element with a 'reference' parameter set for 
aj values above the transition from a regular to a quasichaotic regime. Oscillation power of m^ v (a) and mf (b) magnetization 
components together with the magnetoresistance oscillation power (c) calculated from the relation are presented. A region 
on the [aj — /)-plane corresponding to the spectral peaks shown in Fig. (regular precession regime) is marked as a white 
rectangle in the panel (b) 
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FIG. 3: Magnetization dynamics for a nanoelement with the exchange constant reduced down to A — 2.0 ■ 10~ 6 erg/cm (other 
parameters are as for a 'reference' parameter set) represented in the same way as in Fig. Q Note that the maximal oscillation 
power in a regular precession regime is much smaller than in Fig. 
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FIG. 4: The same as in Fig. 0for a nanoelement with the 'reference' parameter set and the random magnetic anisotropy 
corresponding to a polycrystalline fee Co film structure with the average grain size (D) = 10 nm. 
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FIG. 5: (Color on-line) Comparison of the magnetoresistance power spectra in a quasichaotic regime for nanoelements with 
the 'reference' parameter set (a), the 'reference' parameter set except the reduced exchange stiffness constant A — 2.0 • 1CP 6 
erg/cm (b), the 'reference' parameter set and a polycrystalline structure with average grain size (D) = 10 nm) (c). For the last 
case the random orientations of anisotropy axes for different grains and the fee Co structure (cubic anisotropy of grains with 
K1 uh = 6.0 -10 5 erg/cm 3 ) were assumed. It can be seen that both the reduction of the exchange stiffness constant and inclusion 
of a random anisotropy lead to (i) the expansion of the aj region where quasichaotic oscillations exist and (ii) decrease of the 
out-of-plane oscillation power. 



FIG. 6: (Color on-line) Maps of the magnetoresistance oscillation power for a nanoelement with the 'reference' parameter set 
and a polycrystalline structure with hephep Co grains. The grains with (D) = 10 nm have a uniaxial anisotropy (Ki n — 4.5- 10 6 
erg/cm 3 ) and randomly oriented anisotropy axes. Maps (a)-(d) represent oscillation power spectra for nanoelements with the 
same macroscopic parameters, but different realizations of a random polycrystalline structure. 
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FIG. 7: Magnetization trajectories (a,d) and oscillation power spectra for mj v (b,e) and mf (c,f) magnetization components 
for a random polycrystalline structure realization marked as (b) in Fig. El The graphs in the upper row are for the current value 
aj = 0.50 before the transition from a regular to a quasichaotic behaviour, the lower row shows results for aj = 0.54 which is 
immediately after this transition. Maps in the graphs (b),(c),(e) and (f) display the spatial distribution of the oscillation power 
for the corresponding magnetization components. 
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FIG. 8: (Color on-line) Magnetization dynamics of a polycrystalline (fee Co) nanoelement at T = 300 K: transition from a 
regular to a quasichaotic behaviour (a) and the map of the magnetoresistance oscillation power for the whole aj-range (b). 
Comparison of the latter map to the analogous map for T = in Fig. ^c) shows that thermal fluctuations lead to (i) an 
appearance of an extended region of a regular precession and (ii) a strong power decrease of the 'out-of-plane' oscillations. 
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